STATICS OF POINT JOSEPHSON JUNCTIONS IN A MICRO STRIP 

LINE. 
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Abstract. We model the static behavior of point Josephson junctions in a micro strip line using 
a ID linear differential equation with delta distributed sine non-linearities. We analyze the maximum 
current j max crossing the micro strip for a given magnetic field H. In particular we establish its 
periodicity and analyze how it is affected by the geometry, length, type of current feed, position and 
area of the junctions. For small currents, which is the rule in practice, we show that y m ax can De 
obtained by a simple formula, the magnetic approximation. This model is in excellent agreement 
with measurements obtained for real devices. 
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1. Introduction. The coupling of two low T c superconductors across a thin 
J> ■ oxide layer is described by the Josephson equations pQ. 
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where V and / are, respectively, the voltage and current across the barrier, s is the 
contact surface, J c is the critical current density and $o = 7i/2e is the reduced flux 
quantum. These two Josephson relations together with Maxwell's equations imply 
the modulation of DC current by an external magnetic field in the static regime and 
the conversion of AC current into microwave radiation [21 E]. Other applications 
include Rapid single flux quantum logic electronics and microwave signal mixers 
used in integrated receivers for radio-astronomy 0. In all these systems there is a 
characteristic length which reduces to the Josephson length, Aj, the ratio of the 
electromagnetic flux to the quantum flux $o for standard junctions. 

For many applications and in order to protect the junction, Josephson junctions 
are embedded in a so called microstrip line which is the capacitor made by the overlap 
of the two superconducting layers. This is the so-called "window geometry" where 
the phase difference between the top and bottom layer satisfies an inhomogeneous 2D 
damped driven sine Gordon equation resulting from Maxwell's equations and the 
Josephson constitutive relations 1)1. The damping is due to the normal electrons 
and the driving through the boundary conditions with an external current or magnetic 
field applied to the device. 

Even in the static regime the 2D problem is complicated because of the multi- 
plicity of solutions due to the sine term. However flux penetration occurs along the 
direction of the magnetic field so one direction dominates the other. A quantity mea- 
sured by experimentalists is the maximum (static) current I m & x (H) that can cross 
the device for a given magnetic field H. This gives informations on the quality of the 
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Fig. 1.1. The left panel shows the top view of a superconducting microstrip line containing 
three Josephson junctions, H, I and tf> are respectively the applied magnetic field, current and the 
phase difference between the two superconducting layers. The phase difference <f> between the two 
superconducting layers satisfies —A<f> = in the linear part and — Arp + sin(</>) = in the Josephson 
junctions. The right panel shows the associated 2D domain of size I X w containing n = 3 junctions 
placed at the positions y = w/2 and x = ai, i = 1, n. 



junctions. An important issue is how defects in the coupling will affect this maxi- 
mum current. In particular high T c superconductors can be described as Josephson 
junctions where the critical current density is a rapidly varying function of the po- 
sition, due to grain boundaries. Fehrenbacher et al^3] calculated I ma3L (H) for such 
disordered long Josephson junctions and for a periodic array of defects. Experiments 
were also done by Itzler and Tinkham on large 2D disordered junctions However 
the overall picture is complex and it is difficult from the curve Imax(H) to obtain 
geometric information on the junction. The analysis of such a 2D problem pro- 
vided bounds on the gradient of the solution that were independent of the area of the 
junctions so that little information could be obtained on I max (H). However the study 
[TTj proved the existence of solutions and the convergence of the Picard iteration to 
obtain them. 

Small junctions of length Wi < Xj are easier to study and lead to the well known 
Ima,x(H) = sin(iJwi)/i?[2j. Two such junctions are commonly associated to form a 
Superconducting Quantum Interference Device (SQUID) now routinely used to mea- 
sure magnetic fields. More junctions can be used to form arrays ^1 that can bear 
more critical current and are more flexible than a long junction because the area of 
the junction components and their position can be varied. When the junctions are 
closer than Aj, such arrays behave as a long junction and could be used as microwave 
generators. Almost all models are discrete lumped models where the effect of the 
space between junctions is neglected. In particular the interaction of the junctions 
through this passive region has always been neglected. This makes it difficult to de- 
scribe junctions of different areas, placed non uniformly in the microstrip. This is 
why up to now mostly equidistant and identical junctions have been considered. To 
overcome these difficulties we recently introduced a continuous/discrete model that 
preserves the continuity of the phase and its normal gradient across the junction in- 
terface and where the phase is assumed constant in the junctions. The ID dynamics 
[U] of one junction in a cavity revealed that the junction could stop waves across the 
cavity or enhance them throughout. In |H] this model was used to calculate I max .(H) 
for a miss-aligned SQUID in a 2D cavity. Nonuniform arrays of junctions that are 
generalized SQUIDs have been produced and analyzed in particular by Salez et al 
at the Observatory of Paris 4 and our analysis is in excellent agreement with the 
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measured I max (H). 

In this article we will concentrate on the ID static problem and show that it 
allows for an in depth analysis that was out of reach in the general 2D case. In 
particular we will show the properties of J ma x(-ff), its periodicity, its regularity, the 
relation between different types of current feeds and how it is affected by the posi- 
tion of the array in the microstrip. In addition we introduce and justify the so-called 
magnetic approximation where many details of Imax(-ff) can be controlled. Specifi- 
cally in section two we introduce our model and give preliminary analytical results in 
section three. Section four details the intrinsic properties of the maximal current as 
a function of the magnetic field: its periodicity, the relation between the inline and 
overlap current feed and the simple magnetic approximation. Section five introduces 
two numerical ways to solve the problem. In the sixth section, we study a SQUID 
and examine the effect of a little difference between the junction parameters and we 
compare this to the experiment. Section seven deals with devices with more junctions, 
there we analyze the effect of separating one junction from the others and show the 
agreement with the experimental results. 

2. The model. The device we model shown in Fig. II. H is a so-called microstrip 
cavity (grey area in Fig. 11.1(1 between two superconducting layers. Inside this mi- 
crostrip there are regions where the oxide layer is very thin (~ 10 Angstrom) enabling 
Josephson coupling between the top and bottom superconductors. The dimensions of 
the microstrip are about 100 /im in length and 20 /im in width. The phase difference 
between the top and bottom superconducting layers obeys in the static regime the 
following semilinear elliptic partial differential equation [5] 



(2.1) 



Aip + g(x,y) sin ip = 0, 



where g{x,y) is 1 in the Josephson junctions and outside. This formulation guaran- 
tees the continuity of the normal gradient of tp, the electrical current on the junction 
interface. The unit of space is the Josephson length Aj, the ratio of the flux formed 
with the critical current density and the surface inductance to the flux quantum (f> - 
The boundary conditions representing an external current input / or an applied 
magnetic field H (along the y axis) are 
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= 1 shown in Fig. 11.11 

0, w is called overlap 



where < v < 1 gives the type of current feed. The case v 
where the current is only applied to the long boundaries y - 
feed while v = corresponds to the inline feed. 

We consider long and narrow strips containing a few small junctions of size Wj x un- 
placed on the line y = w/2 and centered on x = a,, i = l,n as shown in Fig. 11.11 
We then search ip in the form 
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where the first term takes care of the y boundary condition. For narrow strips w < 7T, 
only the first transverse mode needs to be taken into account E] because the 
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curvature of ip due to current remains small. Inserting l|2.3[l into (|2.1|) and projecting 
on the zero mode we obtain the following equation for <j) where the O's have been 
dropped for simplicity 
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(2.4) -<j>" + g(x,^ ^sin</.= v- 

where 7 = I/w and the boundary conditions 4>'(0) = H — (1 — v)j/2 > and 4>'(l) = 
H+ (1-^)7/2. 

As the area of the junction is reduced, the total Josephson current is reduced and 
tends to zero. To describe small junctions where the phase variation can be neglected 
but that can carry a significant current, we introduce the following function 

(2.5) 9h( x ) = ^ f° r a-i — h < x < a,i + h, gh{%) = elsewhere, 
where i — 1, . . . n. In the limit h — > we obtain our final delta function model [Sj 

n 

(2.6) — 4>" + djS(x — aj) sin (f) = vj, 

i=l 

where 

(2.7) d, = ^, 3 = j 
and the boundary conditions are 

(2.8) 0'(O)= H -(l-2/) 7 /2, 4>\l) = H + (1 - i/) 7 /2. 

This is our continuous/discrete ID model of a parallel array of many point Josephson 
junctions embedded in micro strip cavity. It preserves the spatial degrees of freedom 
in the linear cavity and the matching conditions at the junction interfaces. 

3. General properties. The delta function seems to be a theoretical way to 
approach the problem. Nevertheless we will show that it provides an excellent agree- 
ment with experiments, in addition to allow simple calculations. We have the following 
properties. 

1. Integrating twice H2.6|l shows that the solution 4> is continuous at the junctions 
x = a,i, i = 1, . . . n. 

2. Let be a solution of the equation (|2.6|l . then <fi + 2kir is also a solution. 

3. Almost everywhere, —cf)"(x) = vy/l, so that outside the junctions, </> is a 
second degree polynomial by parts, 

(3.1) <j)(x) = - ^-x 2 + BiX + d , Vie]fl,,a 1+ i[. 

4. At each junction (x = a^), <f>' is not defined, but choosing t\ > 0, and 62 > 0, 
we get 

lim lim / ' ' (j>"(x)dx = [ * <t>"(x)dx = [(f>'(x)} a t . 
ei^Oe 2 ^oJ a ._ ei J a - °-i 

Since the phase is continuous at the junction x = aj, we obtain: 

(3.2) =*sin(&) , 
with <j)i = <j)(cii). 
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5. Integrating ()2.6)1 over the whole domain, 

r-l n 
W\o = I 4>"dx = y2 d i sin(<^i) - ^7 , 
J o i=l 

and taking into account the boundary conditions, we obtain 

n 

(3.3) 7 = ^ di sin(<£i) , 

8=1 

which indicates the conservation of current. Note that the total current is equal to 
the sum of the jumps of <£>' . 

3.1. Polynomial by part. Let be a solution of 1)2.6)1 and </>i = 4>(<Xi). From 
remark l|3.1|l . tj> is a polynomial by parts. We define Pj+i(x) the second degree polyno- 
mial such that Pi-\-\(x) = 4>{x) for a, < x < dj+i. Using the left boundary condition 
we can specify 4> on [0, oi]: 

(3.4) P 1 (x) = - V ± (x 2 - a\) + (h- {x - a x ) + 0x . 

At the junctions (|3.2I) tells us that Vfc e {1, . . . , n}, 

(3.5) ^fc+iK) - ^-K) = 4 sm(P k (a k )). 

Considering that </>" = — i/j on each interval, the previous relation and the continuity 
of the phase at the junction, we can give a first expression for Pk+i, 

(3.6) Pk+i{x) = — —(x - a k ) 2 + [Pfc(afe) + 4sinP fe (a fc )] (x - a k ) + Pfc(ofc). 

Notice that Pfc+i(x) depends on P k (x), is, j and i7. The parameters v and / are fixed 
by the geometry of the device. So by recurrence we see that </> is entirely determined 
by the values of (f>i, 7 and H. 

From 1)3. 5J) we can obtain another expression for Pfc+i 

(3.7) Pk+i(x) - Pk{x) = d k sm(P k (a k ))(x - a k ). 
Summing all these relations yields 

fc 

(3.8) P k+1 (x)=P 1 (x)+^2d i wa.{P i (a i )){x-a i ). 

i=i 

Polynomials (|3.4(l and 1)3. 6)1 show by construction, that the constants H, j and <fii, 
determine completely the solution of 1)2.6(1 if it exists. In same way, we can show that 
the three other constants, j, <t > '{ a l) an d 4>i nx <f>- From 1)3.8)1 . we give an expression 
of 4> 

n 

(f>(x) = P 1 (x) + ^'H {x > ai }d i sin(0i)(x - a t ), 

2 = 1 

f 1 x > a- 

where T~t{ x > ai } = < „' — " is the Heaviside function. 
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3.2. An n th order transcendental system. Another way to solve l|2.6|) for 
<fi is to write it as a coupled system of n transcendental equations. For that, we first 
eliminate the constant term by introducing such that 

2 

= 0- v-— =tp- f(x) 

and obtain 

n 

(3.9) - 0" + d ^ x ~ a *) sin ^ " /( a *)) = °> 

with the boundary conditions 

0'(O) = H - (1 - i/)7/2, 0'(Z) = H + (1 + i/)7/2. 

To simplify the notation we will write /, = /(ai) and = ip(a,i). Integrating (|3.9|l 
over the intervals [0, a.^], [a^,a^], .. we obtain the relations 





- mf 


+ di sin(V>x - A) = 0, 




a i 


+ d 2 sin(0 2 - A) = 0, 


(3.10) 




+ d 3 sin(V>3 - h) = 0, 






+ d A sin(V>4 - fa) = 0, 




-wit 


+ d 5 sin(0 5 - / 5 ) = 0, 


where we have assumed n 


— 5 as an example. Now we can use the fact that 0' 



the intervals between the junctions and the boundary conditions to obtain the final 
system 

g-(i-„)Z-^-^ 

2 02 — d\ 

03 - 02 -02 - 01 

«3 ^ «2 &2 — O-l 

(3 11) __ " __ | ^3~^2 

a 4 - a 3 a 3 - a 2 

05 ~ 04 04 - 03 

a5 — a4 a4 — 03 
2 0,5 — a 4 

We will use this formulation as well as the one in the previous subsection to establish 
properties of the solutions and solve the problem numerically using Newton's method. 

4. General properties of r y m ax(H) for an n junction array. The general 
problem is 

n 

(4.1) - (j>"{x) + Y d i S ( x ~ a i) sin(0) = vj. 



+ d\ sin(-0i — 


h) 


= o, 


+ d2 sin(0 2 — 


h) 


= o, 


+ d 3 sin (03 - 


fa) 


= 0, 


+ c?4 sin (0^4 — 


fa) 


= 0, 


+ d 5 sin (0 5 - 


■fa) 


= 0. 
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Fig. 4.1. 'ymax(H) curve for an inline current feed, v = (continuous line) and overlap feed 
v = 1 (dotted line) for a three junction unit {1, 5/2, 5/2 + 5/3}, with d\ = d2 = dz = 1. So l\ = 3/2 
and I2 = 5/3. 

with the boundary conditions 

0'(O) = if - (1 - i/) 7 /2, = if + (1 - i/) 7 /2. 

Experimentalists measure the maximum current 7 for a given magnetic field if and 
plot this as a curve "f m ax(H). To compare with real data it is therefore important to 
compute and analyze this quantity. In this section, we give some properties of the 
lmax{H) curve. In the appendix some analytical estimates on the influence of the 
geometry on the maximal current will be presented. 

4.1. Periodicity. We introduce 

lj = aj+i - a,j, 

the distance between two consecutive junctions. Let l m i n be the smallest distance lj. 
We define the array as harmonic if U is a multiple of l m m for ah *• 

Proposition 4.1 (Periodicity of the device) . For a harmonic array, thej max (H) 
curve is periodic with a period 2ir/l m i n . 

Proof. Let <j> be a solution of 14.1fl for a current 7 and a magnetic field H . We 
introduce f(x) = (2tt /l m in){x — a±) and tp(x) = 4>{x) + f(x). So ip verifies 

n 

(4.2) - 4>"(x) + S{x ~ <H) sin(^ - /) = uj. 

i=l 

with V'(0) = H + 2Tr/k - (1 - z/) 7 /2, and ^'(Z) = if + 2tt/^ + (1 - i/) 7 /2. Since, 
/(%) = 2fc7r, Vi G {1, n}, then ^ is a solution of l|4.1|) for if + H p = if + 27r/Z m j n 
and the same 7 , so j ma x(H + H p ) > (if). 

Conversely, by subtracting / from a solution associated to if + H p and a current 
7 , we obtain a solution for if and the same current 7 so j m ax(H + H p ) < j m ax(H). 
From the two inequalities we get 

(4.3) JmaxiH + H p ) = 7mox(if) . 

with ff p = 2Tr/l min . □ 
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In the non harmonic case, if the junctions are set such that lj = Pj/qj, where Pj 
and <7j are integers, prime with each other, then j ma x is periodic with period H p such 
that 

see Fig. 14. II where LCM is the Lowest Common Multiple and HCF the Highest Com- 
mon Factor. To prove this write f(x) = p(x—ax) and use again the previous argument. 
In Fig. 14.11 we show the ~f m ax{H) curve for a three junction unit such that 1% = 3/2 
and h = 5/3 so that the period of 7mox (-f0 is H p = 2irLCM(2,S)/HCF(3,5) = 12vr. 
In the following plots we will only show one period of jmax (H) ■ 

In the general case, we only have an approximate periodicity of 7 „ laa; (H) which 
can be estimated using (|4.4|l . Also, real junctions have a finite size which causes 
lmax{H) — > 0, when H — * +oo. Our model is thus valid as long as the dimensionless 
magnetic field H is not larger than 1/wj. 

4.2. Influence of the position of the junction unit. In this section, we 
examine how the position of the set of junctions in the microstrip (linear domain) 
will affect the j ma x(H) curve. For an array of junctions placed at the distances 
{<Zj, i = 1, n}, we define a junction unit as the set {Zj, i = 1, n — 1}. Then the array 
where the junctions are at {ax + c, a 2 + c, a n + c} is the same junction unit. We 
define a\ as the position of the junction unit. The length of the junction unit is 
lb = a n — a\- The array is centered if (a n + <zi)/2 = 1/2. 

Inline current feed: (v = 0) 

Then the boundary conditions at the edge of the junction unit are <j>'{a^) = 
4>'(0) = H- 7/2, and 4>'{a+) = <p'(l) = H + 7/2, independently of the position of the 
junction unit. 

Proposition 4.2 (Inline junction unit). For inline current feed, jmaxiFI) is 
independent of ax (the position of the junction unit) and of the length I of the cavity. 

Proof. Let 4>x(x) be a solution of Q4.ip . for given 7, H . Let us change the position 
of the junction unit to ai + c so that the junctions are now placed at {ai + c, a 2 + 
c, . . . , a n + c} . It is easy to see that 02 (x) = <px{x — c) satisfies the boundary conditions 
and is a solution. This one to one map between <f>x and <f>2 exists for all c, H and 7 so 
the two junction units have the same 7 mM (-ff). □ 

Then the jmaxiH) curve is independent of the position of junction unit when 
v = 0. By the same argument, we can show that 7 mn i(-ff) is independent of the 
length I of the circuit (see Fig. 14.3(1 . This curve depends only on the junction unit. 

General current feed: (0 < v < 1) 

Then the boundary conditions at the edge of the junction unit are: 

4HflZ) = -vjax +H+(l- i/) 7 /2, (f>'(a+) - H - (1 - i/) 7 /2 + uj(l - a n ). 

Contrary to the inline feed, we cannot shift the phase to find a solution when the 
junction unit has been shifted, because now the boundary conditions depend on the 
position of the junction unit. Consider the derivative $ at the boundaries of the 
junction unit. We will compare the curves j m ax(H) for a centered unit and for a non 
centered unit. For a centered unit, a n — ax = //2 so that 

<t/fc)-H = -(<t/(a+)-H), 
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Fig. 4.2. P/ot o/ "Ymax{H) for a four junctions device (I = 10, di = 1) such li = 1.5, I2 = 2.5, 
£3 = 2. TTte left panel shows a non centered junction unit with a\ = 0.1 and the right panel a 
centered unit with ai = 2. Notice the current dependent shift {4-5f for the overlap solution as one 
goes from a centered junction unit (right panel) to an off-centered junction unit (right panel). The 
junction unit was moved to the left. 



but this equality is false for a non centered unit. It is possible to choose a correction 
H v to the magnetic field H in order to obtain the equality: 

(4.5) H, 

Let us consider two arrays, 1 with a centered junction unit and 2 with the same 
junction unit but non centered. 

Proposition 4.3 (Magnetic shift). Let (H : ~f max ) be the coordinates of a point 
of the j m ax(H) curve for the circuit 1. Then (H + H u ,^ max ) is a point of the curve 
for the circuit 2. 

So, moving a junction unit translates the j m ax curve by vja\. Fig. 14.21 shows 
a lmax{H) for a four junctions device with a non centered junction unit in the left 
panel and a centered junction unit in the right panel. Both inline and overlap current 
feeds are presented. Notice the unchanged behavior for the inline current feed and 
the effect of H u (= — 4.I7/IO) from (|4.5|l in the overlap case. 

Proof. Let </>{jj )7i & 1 } be a solution for an array A\ = {a\, ...,a n } with a centered 
junction unit with 7 and H given. Consider another array Ai with the same junction 
unit moved by s, Ai = {ai + s, a n + s} = A\ + s, the coefficients d\, d n being 
equal for the two circuits. From the solution (^{h^m} ^ or ^ l we can deduce a solution 
^{H+Hv^fii+s} f° r ^2- From (|4.5|l we have: 

^'{H+H„riM+s}( a i + s ) = 9 S '{H,j,6 1 }( a r) 

Taking, 

and from the unicity of the solution, we obtain <f> = ip in the two junction units. Thus, 
if is a solution for {H, 7} given for A%, then ip is a solution for {H + H„, 7} for A2, 
and vice versa. 



H + H„ 



v.) 



- a n ) - (1 - v)- +H U 

l b -I 

— ^ 1" a i 
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Let 7 max i and j ma x,2 be the ^ max curves for the arrays Ai and A 2 . From the 
solutions obtained for A\ we build solutions for A 2 . As a consequence, j ma x,i{H + 
H u ) < jmax^iH). On the other side, from solutions of A 2 we build solutions for Ai, 
then j max ,i{H + H v ) > J ma x,2(H). So, we obtain the equality: 

lmaxs{H + H v ) = ^rnax^iH) . 

□ 

Notice that this equality is independent of the number of junctions. 

4.3. Comparison between inline and overlap current feeds. We now com- 
pare the jmax curves for inline and overlap current feed. For one junction, the prob- 
lem can be solved exactly using polynomials by parts (see remark l3~TJl . We obtain, 
lmax{H) — d\, for all v. For two junctions there is the possibility of d\ ^ d 2 and 
this will change "f m ax{H) qualitatively. Let us study the phase difference between two 
junctions. We use remark f3.il and the boundary conditions to get 

02 - 01 = — 7r(a 2 - ax) 2 + (-P'(ai) + d l sin(0i)) (a 2 - a x ) , 

02-01 a 2 + ax 1-U . 

= -vj \-H —7 + di sm(0i) . 

a 2 — cii I z 

If (a 2 + ai)/2 = 1/2 (the junction unit is also centered), as 7 = jl, (|4.6|l becomes: 

(4.7) *L^ = H -l + d l M<h.). 

a 2 — d\ Z 

Note that we can obtain 14. 7|) from (|4.(i|) with v = 0. We have shown that, 

Proposition 4.4 (Equivalence of all current feeds for a centered SQUID.). For 
a centered two junctions device, all current feeds give the same ^max curve. 

For an inline current feed, v = so that the phase difference 02 —0i is independent 
of the position of the junction unit. This is not true for the overlap feed, where moving 
the junction unit causes a "magnetic shift" as seen above in H v equation l|4.5|l . When 
the number of the junctions n > 3, the ^ ma x curve depends on v. The effect of the 
moving the junction unit on the jmax curve was shown above. So, we can reduce the 
study to a centered junction unit. In this case, we have a\ = (I — h)/2, and 

0>r) = 0'(o) + ^ 1 -v^dx, 

(4.8) =H -{l- v )l-vl+ V ±l = H-(l- V -±\l> 
1 ' V 7 2 2 Z 2 V I J 2' 



^(a+)=if+(l 



vh \ 7 
I ) 2' 



with h — a n — a\. We can write vj — (isli,/l)("//lb), and vlb/l — fJ-- So, equation (|4.1|l 
is equivalent to the system: 

n 

(4.9) - <j>"{x) + **0= ~ sin (0) = f*r > 

with, 

0'(ar) =H-(1- M )7/2, 4>'K) = H + (1 - M ) 7 /2. 
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Fig. 4.3. Plot of 'jmax (H) for the same centered junction unit li =1,22 = 4, ?3 = 3, di = 
<fe = d% = d4 = 1 and different lengths I of the microstrip, from top to bottom I = 8, 16 and 64- 
Notice how the overlap solution tends to the inline solution as one increases I. 



As < v < 1, < < If, /I. Note that It, can be considered as the reference length 
of the device. Also note that if I — > +00 then ft — > and the equation (|4.9|) and 
boundary conditions tend to the situation of inline current feed. Fig. 14.31 illustrates 
this convergence when we increase the microstrip length / for a centered junction unit. 
Notice that the solution for the inline feed is not modified by the variation of length. 
With 1 = 8, we have a maximum difference between the solutions for the overlap and 
inline current feeds. As I increases the solution for overlap current feed tend to the 
solution for inline current feed. We prove this in the appendix 'Convergence of the 
solution for a large length I'. 
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Conclusion: In the appendix, we show that when wy/l tends to the solution 
tends to the one for inline current feed. We can get this by increasing the length I 
or shrinking the junction area, (see in appendix 'Convergence by the coefficient di.'). 
We have three parameters: v, I and "/max- I is determined by the circuit, v comes 
from the 2D model, it depends on the width of circuit. The third parameter, can be 
bounded from above: < 7 < Yli^i- We will see in the next section, what is the 
limit of "f m ax(H) for inline and overlap feeds when di are small. 

4.4. The relation between inline and magnetic approximation. The size 

of the junctions Wi < 1 < w so that di « 1 therefore the jump of the gradient of 

the phase across the junctions can be neglected. This is the magnetic approximation 

where only H fixes the phase gradient. In the previous section, we have shown that 

the solution for inline and overlap current feeds converge to the same jmax (H) curve 

for small di. We will show that this limit is the magnetic approximation. 

a + 

Since [</>'] *_ = di sin(0j) (remark 13.21) and j < J2i di/l, then for small di, then (j> 

tends to the linear function 4>{x) = Hx + c. This magnetic approximation seems crude 
but we show that it approaches the solution for inline feed, see the appendix " Inline - 
magnetic convergence" . There we bound the difference between the 7 curves for the 
inline feed and the magnetic approximation. Fig. 14.41 illustrates this convergence as 
di decreases. 

This approximation gives very good results because we work on very small junc- 
tions and the corresponding di » 10 -2 (compared with the values taken in Fig. 14. 4f) . 

4.5. Magnetic approximation. The magnetic approximation is very interest- 
ing because it gives an analytic expression of jmaxiH) and is independent of the value 
of the current and of the scale of the circuit. Here we consider that <f){x) = Hx + c 
and from (|3.3|) 



Notice that c is the only parameter which can be adjusted to reach the maximum. 
To find the "f m ax{H) curve of the magnetic approximation, we take the derivative 



n 




(4.10) 



5c" 




i=l / \i=l 



The values of c canceling dj/dc are 



(4.11) 




and using |4.10p . we have the solution: 



(4.12) 




This jmax curve is a function of H. A similar expression was given by Miller et al 
|10j for homogeneous arrays. Here we generalize this approach to nonhomogeneous 
arrays and justify it rigorously. 




Fig. 4.4. Plot of ^max(H) for the same junction unit and different coefficients di which are 
all equal, from top to bottom di = 1,0.3 and 0.1. The distances between the junctions are l\ = 1.5, 
l 2 = 2.5, l 3 = 2, I = 10. 



Remark: 



If di = d, Vi 6 {1, n}, we can simplify: 



( SiLi cos(Hai) \ 
W = arcta^ E ^ sin(ffa ^ 
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In the same way, 



'Ymax — 



d cos(c max ) I ^sin(-ffaj) J + sm(c max ) I cos(Hai) 




Here changing the value of d will change linearly the amplitude of ^f m ax curve, this 
is not the case for the solutions of the boundary value problem ((2.61) . We can notice 
too, that the ^maxiH) obtained from this approximation is invariant by the transfor- 
mation, 



We will show in the next sections that when di << 1, (|4.12() provides a good 
estimate of the "f m ax(H) curve of a circuit. In addition, from its invariant properties 
we can compare different models and estimate the parameters of the circuit. It is 
a good approximation for the physical device. A cooperation has begun with the 
LERMA at the Observatoire de Paris to match theory and design for this type of 
circuit with specific properties [15). 

5. Numerical solutions. We used two different methods, a stepping in the 
(H, 7) plane using a Newton iteration and what we call the method of implicit curves 
to find the maximal current of Eq. 1(4.1(1 for H given. 

5.1. Newton's method. We start from the system of nonlinear transcendental 
equations (|3 - 1 1 1> which is written for n — 5. Introducing the vector X = (</>i, <f>2 1 ■■■4>n)i 
(13.11(1 can be written as F(X) — where F is a nonlinear map from R n to R n . To 
solve numerically this equation, we use the Newton method. 



where VF(Xk) is the gradient of F evaluated at X = Xk- A first problem is to 
choose the initial vector Xq. For that consider H = 0, there we expect a solution 
such that 7 ~ Y^i d% consequently w 7r/2 [2tt]. We have our initial vector. After 
finding the solution, we step in H and take as initial guess, the previous solution 
found, which for a small step in magnetic field is assumed to be close to the one we 
are looking for. By this way, we obtain a solution with a magnetic field H + dH and a 
current 7. We can then increase 7 until the method does not converge and this gives 
the maximum current 7 max (i? + dH) for increasing H. Similarly we can compute 
JmaxiH) by starting with a large magnetic field and decrease H to 0. This curve will 
in general be different from the one obtained when increasing H due to hysteresis. 
The two curves need to be overlapped to see where is 7 mo3; (ff). So, we introduce 
another method to be sure to obtain directly ^max curve. 

5.2. Implicit curves method. The polynomials 1(3.4(1 and l(3.fi(l establish the 
existence and value of <j> at the junctions. This function should satisfy the boundary 
conditions. The first one 




X k+1 = X k - {VFiX^FiXk), 



ct>>(0) = P{(0)=H- (l-i/) 7 /2 



is true by construction, the second (for n junction circuit) is: 



(5.1) 



P' n+1 (l) = H+{l-v)± 



3 
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Jt/2 



H 



3ti/2 



2it 



Fig. 5.1. Comparison between the Newton method and the implicit curve method for the -f ma x 
curve for a three junction unit ai =1, 122 = 2, 123 = 3,di = d'2 = d^ = 1, v = 1 and I = 10. The 
squares (resp. the +) symbols correspond to the Newton results for decreasing (resp. increasing) H 
and the continuous line corresponds to the results of the implicit curve method. 



is true only for the solutions of Eq. I|4.1|1 . As we have remarked in the section "Poly- 
nomial by part", <f> is entirely determined by <pi, 7 and H. For H given, the solutions 
of Ea. lj5.lj) define a relation between <f>i and 7. So, the maximal current solution 
depends on <j>\ and 7 and Eq. (j5.1jl is the constraint it should satisfy. As the solutions 
<j> are defined modulo 2n, see Ij3.1j) . we can assume <f>i 6 [— tt, 7r]. On the other hand, 
because of Ij3.3j) . 7 G [0, J^i^i]- To solve this problem with Maple. 15], we plot the 
implicit function (the constraint) of the two variables <f>\ and 7 with H and v fixed, 
defined by 

(5-2) P' n+1 \ x=l {4> 1 , 1 ,v,H)-H-^ 1 = 0, 

with (0i,7) £ [ — 7r, 7r] x [0, J2idi]. The program searches, in an exhaustive way, 
the biggest value of 7 of this implicit curve. Incrementing H, we obtain the relation 
'Jmax(H)- We give an expression of P n+ i for two and three junctions, in the Appendix: 
'Implicit curves'. 

Compared to the Newton method detailed in the previous section, this method 
has the advantage to converge to a global maximum -f max , as long as we give enough 
points to plot the implicit curve. Fig. 15. II compares j ma x{H) using the two methods 
for a three junction unit. The solution given by the implicit curve method is in 
continuous line and superposes exactly with the other two plots. With the Newton 
method we can get trapped in local maxima while the implicit curve method always 
gives the global maximum. On the other hand the Newton method is much faster. 

6. Two junctions. We have seen two methods to solve the problem numerically 
and established general properties. Now let us use these results for an array with a 
few junctions. 

6.1. Same junction strength (d = d\ = cfe). In Fig. 16.11 we plot in the 
left panel ^jmaxiH) of a two junction unit. We find the expected periodicity H p = 
2tt/ (02 — ai), with a maximum for H — in the inline case {y = 0). For the overlap 
feed, we have exactly the inline curve plus a magnetic shift. Notice that for the 
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Fig. 6.1. PZot o/ i/ie ~imax curve for a two junction unit such that l\ = 2. In the left panel, 
d\ = d,2 =0.5 while on the right panel d\ = 0I2 = 3. 



inline feed the amplitude of the j m ax curve is not proportional to di, contrary to the 
magnetic approximation. The larger the di the further away the "f m ax (H) curves are 
from the ones given by the magnetic approximation. This is expected because the 
magnetic approximation neglects the effect of di on the phase. In this section, to 
simplify the discussion we will restrict ourselves to the inline current feed. However 
the results will be valid for the general case. The maximum of ^ m ax corresponding to 
H = mod. Hp is the only case where (02 — 0i)/(fl2 — o-i) = H. On the other hand, 
by construction, in the magnetic approximation (02 — 0i)/(a2 — ai) = H for all H. 
In the general case, the closer H is to 7r/(d2 — ai), the further (02 — 0i)/(a2 — a-i) 
is from H. This can be seen in the right panel of Fig. 16.21 So, there will be more 
tunneling current in one junction than in the other. This phenomenon increases as 
H increases from to it / {a 2 — cii). For that value, we have two possible solutions for 
Imax as shown in the left panel of Fig. 16.21 for H = ir/2. 

As the field crosses n/(a 2 — cti) the two junctions behave in opposite fashion as 
shown by the switch of the jumps in 4> x at the junctions, see rig ht panel of Fig. IO 
These two solutions or reversing behavior of junction, cause a jump in "f' max (H). 
As long as the evolution of 4>\ (or ^2) is continuous there is no jump in j' rnax . To 
summarize, the smaller d is, the closer (02 — 0i)/(ct2 — ai) is to H. Another way of 
relaxing this constraint on (02 — 0i)/(fl2 — ai) for a constant d, is to separate the 
junctions and we can show that lb = ai — o-i - * +00 then jmaxiH) — > d\ + di. 



6.2. Regularity of jmaxiH). Junctions are never perfectly similar, small dif- 
ferences in their areas or their critical currents will affect ^max- In the left panel of 
Fig. 16.31 showing ^maxi.H) for a two junction device there is no discontinuity of the 
slope of the curve r y max (H) labeled "non equal", dj ma x/dH exists everywhere. In 
this case, the value of 4>\{H) and 4>2{H) associated to ^maxiH) vary continuously. 

To show this, consider a circuit such that d\ > d2- Fixing 0i = 7r/2, implies 
that d\ — c?2 < 7 and consequently d\ — di < jmax- From remark i|ri.3|l we have 
Imax < di + d 2 . Combining these inequalities we get d\ — d 2 < jmax < di+d 2 . Now 
let us find the values of H for which these bounds can be reached. From remark (|3.3|l 
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Fig. 6.2. Left panel: plot of the phases <f>i and <f>2 as a function of the magnetic field H for the 
same junction unit as the one shown in the left panel of Fig. 16'. 11 On the right panel, we plot <f>( x ) 
for the same device for H < 7r/(a2 — ai) (continuous line) and H > 7r/(a2 — cti) (dashed line). 
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Fig. 6.3. Plot of "Ymax(H) curves for two two junction units with inline current feed. On the 
left panel we compare ■y m ax{H) given by the model for the cases di = d^ and di ^ cfo. The right 
panel shows the fit of the experimental data from a two junction unit of the Observatory of Paris 
(Reproduced with permission of Faouzi Boussaha and Morvan Salez). 



and equation l|3.6(l we have 

7 = d\ sin(</)i) + d 2 sin(^ 2 ) , 



(6.1) 



"J,2 



1-1/ 



— l x + [H - [vax H — I ) j + di sin cj) t )h + fa , 



with lx = a,2 — a\. By substituting the second equality into the first, and taking the 
derivative with respect to fa, we obtain 



(?7 

— — = di cos fa + d 2 

dfa 



a,2 + a\ 1 — v\ dj 
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— + di cos fa }h + l 



COS 02 , 
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Since we search for the maximum of 7, then dj/dcf>i — 0, so that 
(6.2) d\ cos^i = — di (dih cos(^i) + 1) coscj) 2 . 

When 4>i = 7r/2, this condition gives (f> 2 = n/2 modulo w. Now, inserting these 
solutions in (|6.1|) . we obtain the values of H for which these solutions are possible. 



lmax{H) 


H 


di + d 2 


2kn/h + [v{ax + h/2) + (1 - v)l/2](di + da) /I ~ di 


d\ - d 2 


(2k + l)n/h + [i/(o a + h/2) + (1 - i/)i/2](di ~ d 2 )/i " d x 



This enables from the curve ~/max(H) to estimate di and c^. 

We now proceed to give the condition between di and d 2 such that the behavior 
of the jmax relation changes. Since 4>i(H) varies continuously, cos(4> 2 ) takes all the 
values between —1 and 1. We assume that V0i, dih cos(<^i) + l > where h — a 2 — ax 1 . 
We consider two cases: 

1. cos<^x < 0: as cos</>i > —1 from (|6.2|) we obtain, 

d 2 cos(0 2 ) < - d \ , ■ 
1 — dili 

Since cos(</>2) must take all values between —1 and 1 and d2 > 

(6-3) d 2 < . 

1 — dili 

This is the maximal value than d2 can take compared to di . 

2. cos 0i > 0: as cos^i < 1, for the same reason we obtain: 

(6.4) d 2 > -A— . 

1 + dih 

To summarize d~) max (H) / dH does not vary continuously if 
(6-5) tA^<^< " 



1 - dih ~ ~ 1 + dih 
To illustrate this effect we consider the configuration of a a microstrip with inline 
current feed with two Josephson junctions built by Morvan Salez and Faouzi Boussaha 
at the Observatoire de Paris. The results are shown in Fig. 16.31 The square junctions 
have an area of « l/im 2 , the Josephson length is Aj = 5.6/ito and h = a 2 — ai = 
13/irn (using the junction centers). This gives di = d 2 k 0.0357, h ~ 2.32 if the areas 
are equal. However the experimental data does not go to so that the junctions are 
probably slightly different as expected from l|6.3|) and l|6.4|l . 

0.032969 <d 2 < 0.038923 . 

Only a 10% difference in area is enough to give a regular jmaxiH). From the fit of the 
experimental data (right panel of Fig. I6.3|) we can estimate the areas of the junctions 
as w\ — 0.85255/xm 2 and w 2 — 1.1417/^m 2 . 

As we have seen in the previous section, when the jmax curve does not show 
any spike, it is bounded by di + d 2 and \d\ — d 2 \. From this we can obtain the 
characteristics of the two junctions, their critical current density and area except 
that we do not know which junction corresponds to di and which to d 2 . However if 
the "/max does not have any spikes then we can give the exact area of the junctions 
assuming the critical density current is known. 



1 for small junctions this is not a strong constraint, because since di = w^/w << 1, Wi, w << 1 
and li << 1 are about 10 — 2 . 
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Fig. 7.1. Plot of -fmax{H) curves for a two junction unit ai = 1, a-j 2 (l\ \) and a third 
junction placed at 03 = 5, I2 = 3 (left panel) and 03 = 8, I2 = 6 (right panel). All the junctions 
have the coefficient d{ = 1. 



7. Many junctions. A two junctions circuit is a SQUID and shows a simple 
lmax{H). To obtain specific properties for advanced detectors, experimentalists make 
devices with more junctions. 

7.1. 3 Josephson junctions. When we add a new junction to a circuit with 
two junctions, new oscillations appear on JmaxiH). We cannot predict the amplitude 
of the oscillations but from l\ and l 2 we can estimate the number of oscillations in 
one period i.e. the interval [0, H p ], We intoduce the phase difference for H = 0, 
Aipi = cj>i — (fix. Using Proposition 14.11 we can state that as H goes from to H p , 
4>2 ~ <t>i goes from A(f>2 to A(f>2 + 2irl\/l\ = A02 + 2tt. Similarly 03 — <f>x goes from 
A0 3 to A03 + 2tt(1 2 + h)/h = A</>3 + 2ir(l 2 /h + 1) which becomes A</> 3 + 2n(k + 1) 
if the junctions are placed harmonically so that l 2 = kl\. In that case we expect the 
lmaxi.H) curve to present k + 1 bumps within one period. In Fig. 17.11 the junctions 
are placed in a harmonic way 0,3 — 0,2 = k(a2 — ai), where k — 3 (left panel) and k = 6 
(right panel). As expected we see the 4 intermediate "bumps" in the ^maxiH) curve in 
the left panel and 7 "bumps" in the curve of the right panel. We can see the periodicity 
given by H p = 2n/ (02 — ai) = 2ir/l\, which adds new oscillations. This picture shows 
that the closer the third junction is to the junction unit the fewer oscillations there are. 
Then the oscillations have a larger amplitude. These estimations hold approximately 
in the case of an array with more junctions. 

In other words, when 03 — a\ is large as in the right panel of Fig. 17. H and the left 
panel of Fig l7.2l the shape of jmax curve tends to the one for a two junctions circuit. 
We explain this below. 

7.2. Influence of a faraway single junction for the inline current feed. In 

Fig. 17.21 for each panel, we plot a ^ ma x curve, for a n junction unit, and another with 
the same junction unit plus a far away junction. The ^ m ax curve for n + 1 junctions 
look like n junctions curves to which a shift has been added. Let us evaluate this 
shift. 

Remark that for the junction n, using the notations of Eq. (|3.4J) and (|3.6(l . we 

know that <j) n +i is determined by P n +i- If we increase of e, then </>„+i increase 

of e(a„ + i — a n ). Thus, a variation at (/>„ of e — ±7r/(a„+i — a„) is enough to obtain 
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Fig. 7.2. P/ot of -fmaxiH) showing the influence of a far away junction on a junction unit. 
In continuous line we plot ■y rn , a x (H) for the junction unit only and in dashed line we plot ■y m ax (H) 
for the junction unit together with the far away junction. The left panel shows a two junction unit 
together with a third junction and the right panel a three junction unit together with a four junction. 
The H scale on the bottom of the graphs indicates local minima or maxima of the junction unit 
curves and the scale at the top shows these locations shifted by the quantity |7. the continuous 
line corresponds to the junction unit only For all devices di = 1. 



sin</>„ + i = 1. The farther the last junction, the smaller e, and consequently this 
junction has the smallest action on the junction unit. So, in the search of "f m ax, the 
value of sin </> ra+ i is near 1. The 7^aa: curve of a circuit with n + 1 junctions is close 
to 7™ la2 . + d n+ i i.e. the curve for the n junction circuit with n junctions plus the 
maximal contribution of the last junction. 

Now let us assume that sin n +i = 1. Let us recall the boundary conditions of our 
inline current feed problem: <j)'\ , Q m = H =p 7/2. Therefore the boundary conditions 
at the junction unit are </>'(0) = H — 7/2 and for x such that a n < x < a, i+1 , 
4>'(x) = H + 7/2 — d n+1 . As we have done in the section "magnetic shift", let 
H' = H — d n+ i/2. The previous boundary conditions become 

We find the desired boundary values. Finally we obtain: 



(7.1) 



lim 



v n+l 



->+oo 



H + 



d„ 



+1 



,XH) + d n+ i 



Fig. 17.21 illustrates this convergence. 

This argument can not be extended simply to the overlap or general current feed 
for two reasons. First introducing or taking out the last junction a n +\ induces a 
variation of the magnetic shift H v given by (|4.5II . We could estimate it but we have 
the problem that the curvature of </>, for n junctions device is vj jl where j depends 
on the number of the junctions. This will affect the shift between the junctions and 
consequently the curve 7 max - 

However, numerical simulation show that Eq. 1)7.1(1 remains a good approximation 
for the general case (same order from inline) even with a small number of junctions. 
The general feed and inline feed problems coincide when dn+i/^^di tends to zero. 
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Fig. 7.3. Experimental Ims,x{H) for an array of five junctions in a 2D micro strip line built 
by M. Salez and F. Boussaha of the Observatory of Paris (Reproduced with their permission) . The 
measured data is presented by the + symbols and the magnetic approximation result is in continuous 
line. 

Going back to the physical device, this means that the forces of the junctions are very 
small di « 10~ 2 and for these values the inline and overlap results are practically 
indiscernible from the magnetic approximation. Then (|7.1|l can be used. 

7.3. A real device with 5 Josephson junctions. We have compared our 
theory to the experimental results for a device with two Josephson junctions. The 
same team at the Observatoire de Paris, has made a device with five junctions. Here 
the "f m ax curve obtained is totally different from the one for a simple SQUID. The 
parameters are l\ = 20, I2 = 42, Z3 = 12 and I4 = 6. Fig. 17.31 shows the "/max curve 
where the current and magnetic field have been scaled using approximately the same 
factors as for the SQUID of Fig. Ifi.21 Our modeling approach also gives excellent 
agreement for experimental uniform arrays of 5, 10 and 20 junctions. 

8. Conclusion. We have analyzed mathematically a new continuous/discrete 
model for describing arrays of small Josephson junctions. Compared to standard 
"lumped" approaches, we do not approximate the equations, except for neglecting the 
phase variation in the junction. In particular our approach preserves the matching at 
the interface. 

We establish the periodicity of the j m ax(H) curve, show how it depends on the 
position of the array with respect to the microstrip. This is particularly interesting 
to estimate the proportion of inline current feed versus overlap feed. We show how 
separating a junction from an array will influence ^ max {H). 

We introduce a numerical method for estimating ^ m ax{H) which is more reliable 
than the standard Newton method used up to now. 

The relative simplicity of the model allows in depth analysis that is out of reach 
for the 2D model. In particular we show that solutions for general current feed tend 
to the solutions of inline feed when vj/l — > 0. All models reduce to what we call the 
magnetic approximation for small d{. 

Our global model gives a very good agreement with experimental curves obtained 
for arrays of up to five junctions. The simplicity of the magnetic approximation allows 
to address the Inverse problem of determining features of the array from 7 maa; (-ff). 
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9. Appendix. 

9.1. Implicit curves. In this part, we give an example, of P^ +1 (x) for systems 
with three junctions. We denote: 
sini = sin(0i), 
d = (d 1 sin( ( p 1 )-^+H-(l-iy)^)(a 2 -a 1 ) + ^, 

D _ ^(aj + l-aj) 2 
U — 21 

Then equations i|3.4f> and (|3.(i[) give 

(9.1) P^(x) =-^+d 2 sm(-D 1 + C 1 ) + diBm 1 +H-(l-v)Z 



P[{x) = -^H + d 3 sin [-I> a + {-d 2 sin(D 1 - C x ) - 

+di sini +H-(l-v)?-\ (a 3 - a 2 ) -D l +C 1 
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(9.2) +d 2 frin(-I>i + d) + d x sim +H - (1 - z/)|. 

This example shows that Pf.(x) is C°° in the variables (j,4>i,v, H,x). In particular 
P' n {l) is C°° in the variables (7, <f>i , v, H) . 

9.2. The current feed factor v. analytical estimates. Equation lj4.9|l shows 
that we tend to an inline current feed when I is large. However we should show that 
the jmax curve tends to the one for the inline feed. 

Lemma 9.1 (Solution). For all <pi and H, there exists a 7 such that equation 
\4-l)j has a solution. 

Proof. As we have seen in the section " Implicit curves method" , it is sufficient 
to solve equation (f . 1|) : P r ' +1 (0 = H + (1 — v)j/2, to find a solution. Let us fix a 
value for 0i with v, H, I given. If, 7 < — J2?=i di, then P^ +l {l) < H + (1 — v)j/2. 
Conversely when 7 > di, we obtain P^ +1 (l) > H + (1 — v)-f/2. 

But by construction, P^ +1 {1) is a function continuous in all its variables, in particular 
7. Thus we have at least one value of 7 in [— 5^ i=1 di, X)"=i di], such that P^ + i(l) = 
H + (1 — v)j/2 so that it is a solution for that value of <j>i.U 

We want to study the variation of j(H) versus the current feed v. At this point, 
we do not consider the j m ax curve. Let us fix <j)\. Using the previous property, we 
know that there exists at least one solution of equation (|4.1|> , and particularly almost 
one 7. Without changing <fii or H, we plot all the possible 7 versus v in Fig. 19. II We 
call this curve "f(v) curve. To plot this "/(v) curve, we use the same parameter as in 
Fig. 14.31 with H = 1.3617 (see top panel, we choose this H because there is a big 
difference between the solution for inline and overlap current feeds). We choose for 
4>i the value found with Maple giving the maximum j ma x for the inline feed. Fig. 19. II 
top panel, for v = 0, confirms the 7 maa; value found in Fig. 14.31 But for overlap the 
maximum current we can obtain is near 0. So, there is another value of 4>\ for "f ma x 
of overlap current feed (</>i « 0.252). 

Let us study 7(1/) curve. By definition, 7 = Y^i=i dishi(<pi). Let <f)\ be a value 
such dj/dv exists, then: 

dj , d(j>i 

dv t— i dv 



(9.3) 



97 



< 



d<j>i 



di 



With <f>i — 4>(a,i), and we note in the following = lim^o </>' ' (ai — e) (the left derivative 
of <fi). Let us make some remarks: as 4>± is fixed, 



= 

!/=0 

?0i _ d_ ( (I vl b \ 7) _ (I vlb\ d~1 t hi 



dv 



dv dv { \2 21 J 2 I V2 21 J dv 41 
using (|3.4(l and (|3.6(l we can begin iteration, 

a / 

(dj_i^_icos(^_i) 





( d 1 






, 


dv 










M = _ 


( d 1 








dv 


Vdv- 


h7 


~T + 


ov 



dv 

cos(<fo_i) + 



dv 
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Fig. 9.1. Each panels corresponds to device of the panels of the Fig. \4-3\ We plot the implicit 
curve 7(i / ) curve for H = 1.3617, <j>i = 1.3897. This coordinates give the maximum of the -y m ax(H) 
of inline of the top panels of the Fig. \4-3\ From top to bottom we increase the length of the device 
and we notice the stretching of the 7(1/) curve (with the coefficient found in equation \4 . y( : If, /I). 



with li — dj+i — a,i. This last equation can be written, for i > 3: 

fe-i 



di 



v=0 



We obtain, 



07 \ a, - a! 



dv 



u=0 



k=2 



d^i+i r „07 , dd>i ^ d(j) k 

k=2 
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K\ = k 



Kl = k 



h 

21 ' 



«i 



vlb 
~2T 



a>i+i — o>\ 



I 



k 
21 



Applying absolute values, we obtain: 

c?7 



(9.4) 



dv 



dv 



k=2 



dv 



We do not need to find exact expression of \d<pi+\/dv\, we know that it is a linear 
combination of \dj/dv\ and I7I and so is \dcf>2/dv\. Using l|9.4|) . we can show by 
iteration that \d(f>i/dv\ is a linear combination of \d"f/dv\ and |7|. Applying this last 
remark to inequality (|9.3() . we obtain that there exists two constants, C\ and C2 such 



C 2 | 7 |, 



C\ and C2 are combination of di, h, K\ and K\. For v and di sufficiently small 
\d \ < 1, then 



dj 




dj 


dv 







(9.5) 



d~j 



di 



< 



Co 



l-Ci 



This last equation implies local continuity of the 7 curve as a function of v. As we have 
seen in section " Comparison between inline and overlap" increasing I is equivalent to 
decreasing the range of v (given by fx). Ve, 3L/ I > L \i < e. This shows the 
convergence of jmax (0 < v < 1) curve to inline current feed when I — > +00. 

9.3. Convergence by the junction coefficient di. We want to show that the 
general case convergence to inline case, for small di. We have shown in the axiom 
of the previous appendix that for H and <pi given, we can find almost one solution, 
whatever v. This show that for the same <pi, we can find a general and an inline 
solution. Let us define: 

1. P^(x), Vir (z]a ni a n+ i[ a solution of inline problem (y = 0) of this circuit, 7' 
the maximal current associated at the value <p\. 

2. P£(x), \/x G]a n ,a„+i[ a general solution (y ^ 0, same I and same junction 
unit), 7 s the maximal current associated at the value 4>\. 

' 07 = Pf( aj )-Pf( aj ), 

v Pi = PfM-Pfai)- 

As Pf(ai) = (j>\ = Pl(a t ), we have Pi = 0. We can calculate ai using JO}, 



3. aj and /3,- by: 



oil 



- P9 



Pf(ai)-Pl(ai) = - 



7 1 — <-y« 



21 1 



But 7* and 7 s are positive, so 
(9.6) I 



1 vl b 
ai\<[- + ^- 



i=l 



The aim of this following part is to bound Pi. We proceed by iteration. We recall 
that Ik = ctfc+i — a-k- Using l|3.6(l we estimate Pk+i- 



(9.7) 



A 



fc+i 



2/ 



$ + [d k (sm(P°(a k )) - sm(P l k {a k ))) + a h ] h + Pk- 
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Let us focus on the sine terms, 
sin 



(P°(a k )) - sm(n(a k )) = sra(Jfta fc ) + &) - sin(P*(a fe )), 

= sin(P^(a fe )) [cos(/3 fc ) - 1] + sin(/3 fc ) cos(P^(a fc )) 



compared to ft. From (|9.7() . 



< lft| 2 + lftl 

ain the equivalences sin(ft) ss ft and cos(ft) — 1 p 
ign of sin(P^(a fc )) or cos(P^(a fe )). We neglect |ft| 



We assume ft < < 1 , thus we obtain the equivalences sin( 

-0%, but we cannot predict the sigr " r " '*■•• •*• ••• ' ■ " ' 1 



I/? 



'n+l 



< 



\a n +\\ < 



V l 9 ,2 

21 n 



+ Kilftt + \a n \) in + \Pn\, 

+ d n \(3„\ + \a n \. 



Let us note G = ^X^i=i di/l, we obtain a simple system 
(9.8) Cn+i < M n ( n + GV n , 

Iftil \ „ , r ( d n l n + 1 l„ 



dnln 

\a n \ I \ d„ 1 



£/2 



I, 



with, Cn = I ' "• j • Mr, = ( " '', ,' ) and l",. 

So, we bound |/3„ and \a n \, with ft| and |ai|. 

(9.9) C„ < M n _i(...(M a (MiCi + GV\) + GV 2 )...) + GV n - X 

When d z -> 0, 

1. G ->• then equation iJOJl tend to < M n _i . . . M 2 MxCx- 

2. Af* f J ^ j then, M fc . . . MiMx -> ^ J ^=1 
From the two previous points, we obtain that 



I ft I < N + M(«i-<*i) + o (XJ di J • 

Using (|9~B|l . we have |ai|(ai — <2i) < Ub/2 + f^/ (20) S"=i di, and previous inequality 
become, Vi £ {1, . . . , n} 



(9.10) 



Iftl < |ft| + Oi 



Remember that we have seen at the beginning that ft = 0, H9.10() show that 7 s tend 
to 7*. Since this convergence occurs independently of (fix, we obtain the convergence 
of the "{ m ax curve. 

9.4. Inline - magnetic convergence. We want to show in this part, the con- 
vergence of an inline solution to the magnetic approximation when di « 1. We 
already know that in this case the "f m ax curve for the general current feed and inline 
feed tend to be equal. By this way, we show that for all v, the J m ax curve of Ea. i|4.1[l 
tends to the magnetic approximation when di << 1 . 

We know that the magnetic approximation is given by f(x) = Hx + c max {H). 
Notice that c max does not depend on the value of 7, see Q4.11[l , We are going to 
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compare the magnetic approximation and the inline current feed solution for the 
same geometry. We proceed as in the previous part, we choose 4>i = Ha\ + c max . 
Remember that in the inline case, is a linear function by parts. Vx G]<Zj, ft;+i[, 

P i+ i(x) = (di sm{P l (a l )) + P!(a,i))(x - en) + P l {a i ). 

Let us define: 

on = P!(a i )-f(a i ) = P!(a i )-H, 
pi = Pi( ai )- f( ai ) . 

We obtain that ot\ = —7/2, Pi = and for a n junction circuit a„+i = 7/2. We 
estimate aj+i: 

a t+ i = di sin(P 4 (aj)) + P-(ai) — H = di sh^P^a*)) + on , 
we obtain, an+i = 2}=i dj sm(Pj(aj)) + oti, thus 

n 

(9.11) \a l+1 \<Y,dk- 

fc=i 

We write Pi+i- 

Pi+i = (di sin(P i (a l )) + P' i (a i ))(a t+ i - on) + P,( a ») _ Ha l+ i + b , 
= (di sin(_P 1 (a l )) + P-fa) - H)(a l+1 - a») + Pi(ai) - Hdi + b , 
= a i+ i(a i+ i - di) + f3i . 

Thus if fa = then 

i-1 

(9.12) Pi = 2J Qfc+l(Qfc+l - Qfc) ■ 

fc=l 

Now using the bounds on the a's and bounding the i,'s we get 

(9.13) |A|<nZ 6 ^d fe . 

fe=i 

This shows that the j max of Ea. (|4.1() tends to the magnetic approximation when 
J2k=i dk tends to 0. 



